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Abstract 

The nature of finite temperature transition in QCD is studied on a lattice with 
Wilson fermion. For massless quarks, the transition is smooth for Np = 2, while it 
is of first order for Np = 3 and 6. These results are consistent with the prediction 
based on universality. For the massive degenerate Np = 3 case, clear two states 
signals are observed for m q < 140MeV. From the study of Np = 2 + 1 case, 
namely two massless u-, d-quark and a s-quark with light mass, our results 
suggest that the finite temperature transition in the real world is first order. 

1 Introduction and The Chiral Properties of 
Wilson Action 

The nature of finite temperature transition in QCD is a fundamental 
information for the study of quark-gluon plasma formation in high energy 
heavy ion collisions, and the evolution of the universe at the very early stage. 

In this article we discuss on the point " what is the nature of finite tem- 
perature phase transition with realistic quark masses" . Quarks with much 
heavier masses than the transition temperature will decouple from the dy- 
namics of transition. Because the physical strange quark mass is the same 
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order of magnitude as the temperature of the transition; 100 ~ 200MeV ,we 
should take into account the dynamics of the u-, d- and s-quark properly. 

In this article, we use Wilson fermion action. With this action quarks in 
any number of flavor(iVi?) are treated with a local action. While the chiral 
symmetry is explicitly broken by the Wilson term even for vanishing bare 
quark masses. However chiral properties have been carefully investigated by 
the following axial vector Ward identity [1]. 

2m, < tt|P 5 |0 >= -m v < tt\A 4 \0 > +0(a) (1) 

where A 4 is the fourth component of axial vector current , and P 5 is pseu- 
doscalar respectively. 

From numerical studies it is found that in the confined phase, m q defined 
by Eq.(l) is proportional to m\. It is further found that in the large (3((3 > 5.5 
for Np = 2) region, the values of quark mass are independent of whether it 
is calculated in confined phase or deconfined phase [2], [3], [4]. In the small (3 
region however, m q in the deconfined phase does not agree with that in the 
confined phase and shows peculiar behavior if plotted as a function of 1/K, 
which may be due to 0(a) lattice artifacts [4], [5]. 

Chiral limit (K c (/3)) is defined by extrapolating either m\ or m q to zero 
in the confined region for (3 < 5.5, while for larger (3 region it can also be 
defined directly by the point m q = 0. 

We further define K t} which is the hopping parameter at which the finite 
temperature phase transition occurs. It is a function of temporal size of 
lattice N t} number of flavor Np and (3. As first noticed by Fukugita et a!. [6], 
whether the K t line crosses K c line or not was not trivial. This question 
can be answered if we know whether the point (3 = and K c = 0.25 is in 
the confined phase or deconfined phase when number of flavor is changed. 
The confinement property at (3 = has been studied by Iwasaki et a!. [7]. It 
has been concluded that for Np > 7 the point (3 = and K = K c is in the 
deconfined phase even in zero temperature(we have confirmed deconfinement 
up to N t = 18) and for Np < 6 it is in the confined phase for N t > 4 . 

In this article we present the results of simulations on 8 2 X 10 X 4, 12 3 X 4 
and 12 3 X 6 lattices. In the calculation of hadron mass, the size of third 
direction is duplicated. For the further details of simulations see[8] and [9]. 
In section 2, we show that the nature of chiral transition is continuous for 
Np = 2 and first order for Np = 3 and 6. In section 3 the nature of 
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the transition with massive quarks is discussed for Np = 3. In the case of 
degenerate m g , the first order transition is observed for m q < lAOMeV. For 
Np = 2 + 1 case, namely m u = ~ and m s ^ 0, two states signals are 
observed both for m s ~ lbOMeV and m s ~ AOOMeV. Therefore our results 
indicate first order transition for the real world. The section 5 is devoted to 
conclusions and discussions. 

2 Chiral Phase Transition of Lattice QCD 

In this section, we study the nature of phase transition at m q = 0, namely 
at the point where K c line and K t line crosses. We denote corresponding f3 
as f3 ct . For this study we take a strategy to perform the simulation on K c (f3) 
line starting from the deconfined region and decrease f3 until f3 ct . We call 
this method as "on K c v simulation method. 

It has been shown that in the confined phase, the zero eigenvalues appear 
in the quark matrix around K C} while in the deconfined phase they are not 
small[7]. Therefore the number of iteration for quark matrix inversion Ni nv 
is a good indicator to discriminate the difference of phase on K c . In addi- 
tion to Ni nV} the measurements of Polyakov loops, Wilson loops and hadron 
screening masses are used to discriminate the differences of phases. 

2.1 N F = 2 

Time history of Ni nv for N t = 4 is shown in Fig. 1 . Up to f3 = 4.0 the 
value of Ni nv does not exceed the order of several hundreds, but for f3 < 3.9 it 
increases rapidly with trajectory. Time histories of plaquette, Polyakov loop, 
m\ etc also confirm that the the states for f3 < 3.9 are going into the confined 
phase. Therefore we conclude f3 ct ~ 3.9 — 4.0 for N t = 4. Similar study for 
N t = 6 results in /3 ct ~ 4.0 - 4.2. The K t (/3) line for N t = 4 and K c (/3) 
lines are shown in Fig. 2[6],[3],[10],[11]. The crossing points f3 ct obtained by 
"on K c v simulation method are consistent with the points determined by the 
extrapolation of the K t line toward K C} both for Nj = 4 and 6. 

The m\ on K c in the deconfined phase, are shown in Fig 3. They are 
decreasing gradually and at the region of f3 ct , they are consistent with zero. 
Around f3 ct we could not observe two states signal. From these results we 
conclude that the phase transition in this case is continuous. 



3 



2.2 N F = 3,6 

The same method as in the case of Np = 2 is applied for the study of 
Np = 3 and N t = 4. We observe clear two states signal at f3 = 3.0; if a 
initial condition is deconfined states, Ni nv is stable around several hundreds 
and system stays in the deconfined phase, while if we start from a mixed 
states, Ni nv increases with trajectory. Similarly the hadron screening mass, 
quark mass, plaquette and Polyakov loop also suggest two states signal at 
P = 3.0[8]. 

The situation for Np = 6 is completely same as the case of Np = 3 except 
that two-states signal is observed at f3 = 0.3. 

From these results we concludes that the nature of chiral transition for 
Np = 2 is continuous while for Np = 3 and 6 , first order. These results are 
consistent with those from staggerd quark [12], and the prediction based on 
universality arguments [13]. 

3 Finite Temperature Phase Transition with 
Massive Quark 

3.1 Degenerate Np = 3 Case 

Before proceeding to the realistic non-degenerate quark mass case, we 
study degenerate Np = 3 case. In this case the simulation is made on K t 
line for Np = 3 in Fig. 2. We observe clear two states signal at f3 = 4, 0, 4, 5 
and 4.7(Fig. 4). However no two states signal has been observed at f3 = 5.0 
and 5.5. At f3 = 5.0 we have made more fine tuning of K than in the case of 
f3 = 4.7[9]. 

Then the critical quark mass m " 1 up to which the clear two states signal 
is observed is bounded from below by the m q at f3 = 4.7 and K = 0.1795. 
That is m="*a > 0.175. 

The lattice spacing around this f3 region is determined by the p meson 
spectroscopy. In Fig. 5 the m p a and m q a are plotted as a function of 1/ K — 
1/K C . They are almost independent of Np, N t and f3. From m p a at K c and 
physical p meson mass(m p = 770AfeV), we get a -1 ~ 0.8GeV. Then m c q Tli 
in physical unit becomes > 140AfeV. 
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3.2 N F = 2 + 1 



We proceed to the study of realistic quark mass case: m u = m<i ~ and 
m s ^ 0. As the m " 1 for degenerate Np = 3 case discussed in the previous 
section is almost same order of magnitude as the physical strange quark mass, 
the nature of transition will depends on m s sensitively. 

In this case, confined and deconfined phases are separated by the finite 
temperature transition line of strange quark K s p } which is shown schemat- 
ically in Fig. 6. Along K s p line, the nature of transition changes from first 
order at f3 ~ 3.0 to continuous when f3 is increased. We want to know 
whether the critical quark mass of s-quark m c s rii at which the nature of tran- 
sition changes from first order to continuous is lighter than physical strange 
quark mass or not. 

For this study we use the similar method as "on K c v simulation method; 
namely keeping K u = Kd = K c , and K s to the value corresponding to of 
m s ~ lbOMeV or AOOMeV , and decrease f3 in the deconfining region, until 
we hit K s t line. The K s for m s ~ lbOMeV and AOOMeV are obtained from 
Fig. 5 [9] and schematically shown in Fig. 6. 

In Fig. 7 we show time history of plaquette at the crossing point for 
m s ~ AOOMeV. Clear two states signal is observed at f3 = 3.9, both on 
8 2 X 10 X 4 and 12 3 X 4 lattices. The mesurements of plauette, Polyakov loop 
and m\ etc also indicate two states signal[9]. For m s ~ lbOMeV too, two 
states signal is observed at f3 = 3.5. 

4 Conclusions and Discussions 

Our results implies that chiral transition is continuous for Np = 2 (both 
for N t = 4 and 6), and for Np = 3 and 6 it is first order. These results 
are consistent with those of staggerd quark[12], and the prediction based 
on universality [13]. For massive degenerate Np = 3 case, our result imply 
m c q rtt > 140MeV ,and for N F = 2 + 1 > 400MeV. 

These critical quark masses are much heavier than those of staggerd 
quark[14],[15]. Columbia group reported that first order transition is not 
observed at m u a = m^a = 0.025, m s a ~ 0.1. Using a -1 ~ O.bGeV for 
staggerd quark around this f3 region[9],[16] ,the quark masses in the physical 
unit are m u = m<i ~ 12MeV and m s ~ bOMeV respectively. One of the 
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pssibilities for the discrepancy will be that f3 in the both simulations are 
far from continuum limit. To conclude the nature of transition in the real 
world, further studies with larger f3 are necessary, but our result suggests a 
possibility for the first order transition in QCD in the real world. 

Simulations are performed with HITAC S820/80 at KEK, and with VPP- 
500/30 and QCDPAX at University of Tsukuba. We thank members of 
KEK for their hospitality and strong support. We also thank other mai- 
mers of QCDPAX collaboration for their support. This work is in part 
supported by the Grant-in- Aid of Ministry of Education, Science and Cul- 
ture(No.06NP0601). 
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= 13cm Figl.ps 

Figure 1: Time history of Ni nv on K c for Np = 2 on 8 2 X 10 X 4 lattice. 

= 13cm Fig2.ps 

Figure 2: Phase diagram for Np = 2 and 3. Filled symbols are for A^ c (m 2 ) 
and K c (m q ). Open symbols are for K t (N t = 4) for Np = 2 while open circles 
with cross for Np = 3. Circles are our data. Lines are for guiding eyes. 

= 13cm Fig3.ps 

Figure 3: m 2 a 2 on the K c line for N t = 4 and 6 lattice. The spacial sizes 
are 8 2 X 20 and 12 2 X 24 respectively, where original lattice are duplicated in 
the 3-rd direction. 

= 10. 3cm Fig4.ps 

Figure 4: Time history of the plaquette for Np = 3 at f3 = 4.7 on a 12 3 X 4 
lattice. 

= 12cm Fig5.ps 

Figure 5: m p a and 2m q a in the confined phase. Open symbols are for Np = 2, 
f3 = 3.0, 3.5, 4.0, 4.3, and 4.5 on an 8 2 x 10 x 4 lattice. Filled symbols are 
for N F = 3, f3 = 4.0, 4.5 and 4.7 on 8 2 x 10 x 4 and 12 3 x 4 lattices. K c for 
Np = 2 is used. 

= 12cm Fig6.ps 

Figure 6: Phase diagram for Np = 2+1 as a function of K s at K u = = K c . 
for N F = 2 is ~ 4.0 and for N F = 3 is ~ 3.0 

= 12cm Fig7.ps 

Figure 7: Time history of the plaquette for m s ~ 400 MeV on a 12 3 X 4 
lattice. 
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